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Abstract. Following the classical approach of Polya-Schur theory |14| we 
initiate in this paper the study of linear operators acting on M.[x] and preserving 
either the set of positive univariate polynomials or similar sets of non-negative 
and elliptic polynomials. 
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■ 1. Introduction and main results 

0^ , Let M.[x] denote the ring of univariate polynomials with real coefficients and 

1^ ' denote by K„ [x] its linear subspace consisting of all polynomials of degree less than 

or equal to n. 

In what follows we will discuss the following five important types of univariate 
^ — ^ , polynomials: 

OO i Definition 1.1. A polynomial p(a;) e R[a;] is called 

- hyperbolic, if all its roots are real; 

- elliptic, if it does not have reals roots; 

- positive, if p{x) > for all cc G M; 

- non-negative, if p{x) > for all x S M; 
I - a sum of squares, if there is a positive integer k and there are polynomials 

Pi{x), . . . ,Pk{x) 6 M\x\ such that p{x) = p\{x) + . . . +p\{x). 

Note that the term "elliptic" is sometimes used to define other types of polyno- 
mials, see, e.g., [TU1[T2]. The set of non-negative polynomials is classically compared 
with the set of sums of squares which is a subset of the latter. Moreover, a well- 
known result claims that in the univariate case these two classes coincide, see, e.g., 
[13 p. 132]. 

Proposition 1.2. A polynomial p{x) G M\x\ is non-negative if and only if there 
exist pi(x),p2{x) G R[a;] such that p{x) = p\{x) -['p\{x). 
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Remark 1.3. Note that the situation is quite different for polynomials in several 
variables. In particular, even in 2 variables not all non-negative polynomials can 
be represented as sums of squares. One of the simplest examples of this kind is 
the polynomial p{x, y) — x'^y^ix^ + — 3) + 1 which is non-negative but can not 
be represented as the sum of squares, see [llj for details. In general, this topic is 
related to the Hilbert 17-th problem, see [T^. 

Definition 1.4. Let V denote cither Mjja;] or M\x\. We say that a map ^ -.V 
preserves a certain set M C V if for any p{x) G M its image ^{p{x)) belongs to M . 

In this paper we study linear operators on M.[x\ or Mn[a::] which preserve one of the 
classes of polynomials introduced above. Namely, we call a linear operator acting 
on M\x\ or a hyperbolicity-, ellipticity-, positivity-, non-negativity-preserver 

if it preserves the sets of hyperbolic, elliptic, positive, non-negative polynomials 
respectively. The classical case of (linear) hyperbolicity-preservers which are diag- 
onal in the monomial basis of R[x] was thoroughly studied about a century ago by 
Polya and Schur [T3] . Its substantial generalizations both in the univariate and the 
multivariate cases can be found in [H [51 [3] . 

Following the set-up of [H] we concentrate below on the remaining three classes of 
preservers (restricting our attention mainly to linear ordinary differential operators 
of finite order, see Remark II. 6p . In short, it turns out that there are much fewer 
such linear operators than those preserving hyperbolicity. More precisely, our two 
main results are as follows. 

Theorem A. Let Uq : M[x] M.[x] be a linear ordinary differential operator of 
order k > I with polynomial coefficients Q = (qo{x), qi{x), . . . , qk{x)), qi{x) £ M[a;], 
« = 0, . . . ,fc, qk{x) ^ 0, i.e., 

d d'' 
Uq = qa{x) + qi{x)— + q2{x)-r^ + ... + qk{x)—r. (1) 
ax dx^ dx'^ 

Then for any coefficient sequence Q the operator Uq does not preserve the set of 
non-negative (resp., positive or elliptic) polynomials of degree 2k. 

Corollary. There are no linear ordinary differential operators of positive finite 
order which preserve the set of non-negative (resp., positive or elliptic) polynomials 
in R[x]. 

Remark 1.5. Notice that by contrast with the above situation there are many 
hyperbolicity-preservers which are finite order linear differential operators with 
polynomial coefficients. In fact, such examples exist even among operators with 
constant coefficients, see Rcmark ll.6l and [2]. 

Remark 1.6. Any linear operator on R[x] and C[x] can be represented as a linear 
ordinary differential operator of, in general, infinite order, i.e., as a formal power 
series in ^ with polynomial coefficients, see, e.g., Thus the subclass of finite 
order finear differential operators, i.e., those belonging to the Weyl algebra Ai is 
a natural object of study. Note that unlike the case of finite order operators there 
exist plenty of linear differential operators of infinite order which preserve positivity. 
Apparently, the simplest example of this kind is 

(l- —Y^ - 1 + — — + (2) 

\ dx J dx dx'^ 

More generally, the inverse of any finite order differential operator with constant 
coefficients and positive constant term whose symbol is a hyperbolic polynomial 
yields an example of such an operator. 
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Luckily the case of infinite order linear differential operators with constant co- 
efficients can be handled completely. Namely, slightly generalizing a one hundred 
years old result of Remak [IS and Hurwitz (see also Problem 38 in [15, Ch. 7]) 
one obtains the following statement. 



Theorem B. Let a = (ao, ai, . . . , afc, ) he an infinite sequence of real numbers 
Consider the infinite order linear ordinary differential operator 

d d^ d^ 

' dx^ dx^ 

with constant coefficients. Then the operator Ua preserves positivity (resp., non 
negativity) if and only if one of the following two equivalent conditions holds: 



^ dx 



(3) 



(1) for any positive (resp., non-negative) polynomial p{x) 
one has that 



OkX'^ 



C/q(p(0)) = floOfo + oiai + . . . + klokcek > {resp., > 0) 
(2) the following infinite Hankel matrix 

( ao \-\oL\ 2\oL2 ■ . ■ Uai 



V.ai 
2\a2 



2\a2 
3!q!3 



3!a3 
4!a4 



\ 



{l + l)\ai+i 
{l + 2)\ai+2 

{2l)\a2i 



■J 



is positive definite (resp., positive semi- definite), i.e. all its principal mi- 
nors are positive (resp. non-negative). 

To illustrate the latter result notice that for the operator ^ above one has 
1 = ao = ai = a2 = as = ... and A/ = nj=i(j0^i ^ = 0,1,... where A/ is the 
corresponding (Z + 1) x + 1) principal minor, see 

Remark 1.7. The major remaining challenge in this area is to classify a/Zpositivity- 
preservers. We finish our introduction with this question. 

Problem 1. Find a complete classification of positivity- preservers. 

We also state a more concrete and (hopefully) simpler question. 

Problem 2. Is it true that any positivity-preserver which is an infinite order linear 
differential operator with constant coefficients has a hyperbolicity-preserver as its 
inverse? 

Acknowledgments. The authors are grateful to Petter Branden for important 
references. The second author is sincerely grateful to the Swedish Royal Academy 
of Sciences and the Mittag-Lefiler Institute for supporting his visit to Stockholm 
in Spring 2007 when a substantial part of this project was carried out. The first 
and third authors would like to thank the American Institute of Mathematics for 
its hospitality in May 2007. 

2. Some preliminaries on the considered classes of preservers 

Below we discuss the relationships between the classes of ellipticity-, positivity-, 
and non-negativity-preservers. As we mentioned in the introduction the set of all 
univariate non-negative polynomials coincides with the set of sums of squares and 
therefore linear preservers of the latter set do not require separate consideration. 
On the other hand, it is obvious that the sets of elliptic, positive, and non-negative 
polynomials are distinct. In this section we answer the question about how different 
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are the corresponding sets of ellipticity-, positivity- and non- negativity-preservers, 
respectively, see Theorenis l2.4l and l2.5l below . 

We start with the following lemma showing that the assumption that a linear 
operator $ is a non-negativity-preserver is quite strong. 

Lemma 2.1. Let $ : M.[x] — *■ M.[x] be a linear operator preserving the set of non- 
negative polynomials. //$(!) = then $ = 0. 

Proof. Assume that $(1) = 0. First we show that for any polynomial p{x) = 
dnX^ + . ■ . of even degree n = 2m one has that if a„ > then <i>(p(x)) is non- 
negative and if a„ < then $(p(a;)) is non-positive. Indeed, if n is even and 
a„ > then p{x) has a global minimum, say M. Thus p{x) -I- |Af| > for all 
a; S M. Therefore, ^{p{x) + \M\) > for all a; e K.. However, by linearity and the 
assumption $(1) = we get that 

$(p(a;)) = $(p(a;)) + |M|$(1) = + |M|) > 

for all a; e R. For a„ < the result follows by linearity. 

Now let us show that $(1) = implies that $ = 0. Let q{x) be a polynomial 
such that ^{q{x)) = Onx" + . . . + Oix'^ with the smallest possible non- negative 
value of i such that Oi ^ 0. Let p{x) be a monic real polynomial of even degree 
satisfying the condition (leg{q{x)) < deg{p{x)). Thus p{x) + fJ-q{x) is monic for 
any G R. The above argument shows that ^{p(x) + fiq{x)) is non-negative for 
all /i G M. Notice that our choice of q{x) implies that the polynomial ^{p{x)) has 
vanishing coefficients at the degrees 0, . . . , i — 1. Hence $(p(a;)) — bix^ -f . . . -I- biX^ 
for some positive integer I and some coefficients bi, . . . ,bi G M. Then for any given 
/i there exists (7^(a;) S R[x] such that ^{p{x)) + iJ,^{q{x)) — x'g^(a;). Obviously, 
the constant term of g^{x) equals bi -\- fj,ai. Since ^ there exists fiQ £ M 
such that bi + ^oOi = 3^^0(0) < and by continuity it follows that there exists a 
neighborhood A^(0) of the origin such that g^o{x) < for all x e N{0). Therefore, 
there exists ^ xq £ N{0) such that Xq > 0, hence XQgij,g{xo) < 0. This contradicts 
the assumption that $ is a non-negativity-preserver. Thus ^{q{x)) has a vanishing 
term of degree i, which contradicts the choice of q{x). We deduce that ^{q{x)) = 
for aU e R[a;]. □ 

Theorem 2.2. Let $ : R[x] M.[x] be a linear operator. Then the following 
conditions are equivalent: 

(1) $ preserves the set of elliptic polynomials; 

(2) either $ or — $ preserves the set of positive polynomials. 

Also each of these conditions implies that 

(3) either <J> or — $ preserves the set of non-negative polynomials. 

Proof. Note that the identically zero operator satisfies neither condition (1) nor 
condition (2). Therefore, we will assume that 4> ^ 0. 

(1) ^ (2). Assume that <1> preserves the set of elliptic polynomials and that 
neither <I> nor — $ preserves positivity. In other words, since <& is an ellipticity- 
preserver this means that there exist positive polynomials p{x),q{x) £ M.[x] such 
that ^{p{x)) > and < for all x S R. Note that no elliptic polynomials 

can be annihilated by $ since is not an elliptic polynomial. We consider the 
following two subcases: 

A. There exist two positive polynomials p{x), q{x) as above such that deg ^(p{x)) ^ 
degi>(g(a;)). Wlog we can assume that deg$(p(x)) > deg$(g(a;)). Since ^{p{xj) 
is a positive polynomial it has even degree and positive leading coefficient. Thus 
for any /i £ R the polynomial ^{p{x)) -\- /i$(g(a;)) has the same properties, i.e., is 
of even degree and has positive leading coefficient. Hence there exists xo{^) £ R 
such that -I- > for all x with \x\ > xo{fi). 
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Now set yo := <i>(p)(0) and zq ■= <i>(q)(0). Obviously, yo > and zq < since 
is positive and ^{q) is negative. Let fiQ — > 0. Then p{x) + fJ-Qq{x) is 
positive since it is the sum of two positive polynomials. At the same time for its 
image we have that $(p + fioq)(x) > for x > a;o(/^o)- However, at the origin one 
has 

^{p + Mog)(0) ^yo + ^^zo 2/0 - 2yo = -yo < 0, 
-Zq 

SO by continuity <i>(p + fiQq){x) must have at least one real zero, which is a contra- 
diction. 

B. It remains to consider the case when the images of all positive polynomials have 
the same degree, say m. Let $(p(x)) = amX™' + . . . , ^{q{x)) = bmX™' + . . .. Since 
^{p{x)) > it follows that a™ > 0, and since < one has bm < 0. Thus 

the polynomial —hmp{x) + amq{x) is positive. However, its image is of the degree 
less than m, which is a contradiction. 

(2) =^ (1). If $ is a positivity-preserver then by linearity <& is also a negativity- 
preserver, and thus $ preserves the set of elliptic polynomials as well. 

(2) =J> (3). Assume that $ preserves positivity. Take p{x) e R[a;], p{x) > 0. 
Then for any e > 0, p{x)+e > 0. Thus $(p(a;)) -f = ^{p{x)+e) > 0. Taking 

the limit when e — > we get that > 0. □ 

The following example shows that, in general, (3) does not imply (1) and (2). 

Example 2.3. Let $ ; R[x] R[x] be defined as follows: $(1) = a;^ <I>(x') = 
for all i > 0. Obviously, $ preserves the set of non-negative polynomials but does 
not preserve the set of positive polynomials since 1 is mapped to x^ which is only 
non-negative. 

We are now going to show that in fact this example is in some sense the only 
possibility, i.e., it essentially describes the whole distinction between positivity- and 
non-negativity- preservers. 

Theorem 2.4. Let <I> : M.[x] —^ M[x] be a non-negativity-preserver. Then either $ 
is a positivity-preserver (and therefore an ellipticity-preserver as well) or <&(!) is a 
polynomial which is only non-negative but not positive. Moreover, in the latter case 
for any positive polynomial p{x) G M.[x] the zero locus of ^{p{x)) is a subset of the 
zero locus o/ $(1). 

Proof. Assume that $ is a non- negativity-preserver. Then $ sends positive poly- 
nomials to non-negative ones. Let us assume that p{x) G M[x] is positive but its 
image ^{p{x)) has real zeros. Since p{x) is positive there exists e > such that 
p{x) — e > for all x. Thus for its image we have 

- e) = ^{p{x)) - > 0. 

Set g{x) := $(p(a;)), f{x) := ^{p{x) - e), and h{x) :^ $(1). Since all three 
polynomials are non-negative and g{x) = eh{x) -\- f{x), it follows that for any xq 
such that g{xo) = one has that /(xq) = h{xo) = 0. Since h{x) is a polynomial 
then either h{x) = 0, or h{x) has a finite number of zeros. However, the first 
possibility is ruled out by Lemma since h{x) = $(1). The second possibility 
implies that all positive polynomials whose images are non-negative but not positive 
have altogether only a finite number of zeros belonging to the zero locus of $(1). □ 

Corollary 2.5. Let $ : R[a;] ^ M.[x] be a linear operator such that $(1) > 0. Then 
the conditions (1), (2) and (3) of Theorem \2.Si are equivalent. 

In exactly the same way we can show the following. 
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Theorem 2.6. Let $ : R„[a:] M„[x] be a linear operator with $(1) > 0. Then 
the following conditions are equivalent: 

(1) $ preserves the set of elliptic polynomials of degree < n; 

(2) either $ or — $ preserves the set of positive polynomials of degree < n; 

(3) either $ or — $ preserves the set of non-negative polynomials of degree < n. 

Remark 2.7. Corollary 12.51 and Theorem 12.61 will allow us to reduce the investi- 
gation of non-negativity-, positivity-, and ellipticity-preservers (both in the finite- 
dimensional and the infinite-dimensional cases) to just one of these three classes of 
preservers. 

For the sake of completeness notice that for non-linear operators the situation is 
different from the one above as the following simple examples show. 

Example 2.8. 1. The bijective map $1 : M.[x] R[a;] defined by 

$i(p(a;)) =p{x)+c 

where c is a positive constant, preserves both positivity and non-negativity but does 
not preserve ellipticity. 

2. The bijective map $2 • — > M[x] defined by 



The aim of this section is twofold. Firstly we want to recall what was previously 
known about positivity- and non- negativity-preservers in the classical case of linear 
operators acting diagonally in the standard monomial basis of M.[x] and secondly 
we want to point out some (known to the specialists in the field, [4j, [5]) mistakes 
in the important treatise [8j. 

3.1. Known correct results. Let T^o :M.[x] R[a;] be a linear operator defined 



Denote them by {A^j^Q and {Ai}"^Q, respectively. We will refer to such opera- 
tors as diagonal transformations or diagonal sequences. Diagonal transformations 
preserving the set of positive polynomials are referred to as A-sequences in the lit- 
erature , see [31 E] . Reserving the symbol $ for general linear operators we use in 
this section the notation T G A to emphasize that T is a diagonal transformation 
preserving positivity. Multiplying if necessary all elements of our sequence with 
— 1, we can assume that Aq > 0. 




3. The case of diagonal transformations 




(4) 



(5) 
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Remark 3.1. Notice that in the finite-dimensional case we only need to consider 
transformations acting on M„[a::] for n even since there are no positive polynomials 
of odd degree and a sequence {A^j^^Q^ preserves the set of positive polynomials in 



2k+i if and only if {A^j^^Q preserves the set of positive polynomials in . 



Let us establish some immediate consequences of the fact that a diagonal trans- 
formation T is a positivity-preserver. 

Lemma 3.2. Assume that a transformation Ta = {A^jf^Q, a G NU{cx)}, preserves 
positivity. Then 

(1) Ai > for any even i; 

(2) Xf < AoA2i for any i. 

Proof. To settle (1) consider the polynomial p{x) = -\- 1 which is positive if i is 
even. Thus T{x) = XiX^ + Aq should be positive as well. Since Aq > 0, the result 
follows. 

To settle (2) consider the polynomial p{x) = x^* -|- ax* -|- b with < 46. Then 
p{x) is positive as well as its image q{x) := ^{p{x)) = X2iX^^ +aXiX^ + bXo. If i is odd 
then the positivity of q{x) is equivalent to the negativity of its discriminant, i.e., 
Dq := a^A? — 46AoA2i < 0, which imphes A^ < AoA2i since < 4b. Finally, if i is 
even then q{x) is positive iff either < or > and if additionally both roots 
of the quadratic polynomial A2i2:^ -I- aXiZ + bXo are negative. In the first subcase 
one has A| < AoA2i as for i odd. In the second subcase we obtain that the positive 
polynomial — ax* -I- 6 is transformed to the polynomial A2iX^* — aA^x* + bXo which 
has some real roots. To check this notice that the roots of A2iZ^ aA^z + bXo are 
opposite to that of A2i2;^ -|- aXiZ + bXo and are, therefore, positive. Thus extracting 
their j-th root one will get some positive roots as well. This contradiction shows 
that the inequality A| < AoA2i is necessary for the positivity-preservation. □ 

Diagonal transformations which are positivity-preservers are known to be very 
closely related to the following class of sequences of real numbers. 

Definition 3.3. A sequence {Ai}f^Q is called positive definite if for any positive 
polynomial p{x) = x" + an-ix"~^ -I- ... -I- aix -|- ao G Ma[.x] one has that A„ -I- 
ttn-iXn-i + . . . + aiAi + aoAo > 0, i.e., Ta{p){l) > 0. 

In the infinite-dimensional case the following characterizations of the set of di- 
agonal positivity-preservers is known. 

Theorem 3.4. ([4], O Theorem 1.7]) Let {Xk}kLo be a sequence of real numbers. 
Then the following conditions are equivalent: 

(1) {Xk}f^, e A; 

(2) {Afcj^Q is a positive definite sequence; 
Ao Ai ... Afc 
Ai A2 ... Afc+i 



(3) |(A.+,)| 



> 0, for fc 0, 1,2. 



Afe Afc+i . . . A2fc 
(4) There exists a non- decreasing function fj,{t) with infinitely many points of 
increase such that for all k — 0,1,2, . . . one has 

00 

Xk= t''d^i{t). 



Proof. The equivalences (2) <^ (3) <^ (4) are settled in [T71 p. 132] independently of 
condition (1). The implication (1) (2) is evident so we only have to concentrate 
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on the remaining implication (2) (1). Take a positive polynomial g{x) — aix^ + 

; 

. . . + aix + ao and set q{x) :— A{g{x)) = ^ Xiaix\ We want to show that q{x) > 

1=0 

for all real x. 

oc 

Condition (4) implies that Xi — J f-dait), where ait) is a monotone non- 

— C30 

decreasing function with infinitely many points of increase. Hence 

^ oo OQ ^ OO 

q(x)=J2{[ t'da{t)^x' = I "^a^fx'da{t)= ( g{xt)da{t) > Q 

^=0 - oo - oo '=0 

since g{xt) > for all t. Notice that the above integrals are convergent for any 
fixed value of x. Thus q{x) > and the lemma follows. □ 

In the finite-dimensional case one has a similar statement. 

Theorem 3.5. ^ Let {A^j^^Q be a sequence of real numbers. Then the following 
conditions are equivalent: 

(1) {XkVLo e A; 

(2) {AfcjJJ^Q is a positive definite sequence; 
Xq Ai ... Afc 
Ai A2 ... Afc+i 



(3) |(A.+,)| 



> 0, for fc = 0,1,2. 



Afe Afe+i . . . A2fe 
(4) There exists a non- decreasing function with at least n points of increase 
such that 

00 

Afe = / t''dfi{t), /c = 0,l,2,...,n; 



Proof. Repeats that of Theorem 13.41 □ 

3.2. Known wrong results. To present some erroneous results from 8^ and the 
corresponding counterexamples we need to introduce the following classes of diag- 
onal transformations. 

Definition 3.6. We say that Ta, a G NUoo, or, equivalently, the sequence {Xi}f^Q, 
is a hyperbolicity-preserver, if for any hyperbolic p{x) G Mc([a;] its image Ta{p{x)) is 
hyperbolic. We denote this class of transformations by Ha or H. 

Clearly, this class is the restriction of the earlier defined class of hyperbolicity- 
preservers to diagonal transformations. 

Theorem 4.6.14 of [51 states that T G A if and only if T^^ G H. We wiU now 
show that this statement is wrong in both directions. 

Proposition 3.7. There exist 

(i) T G A such that T-^ H; 

(ii) T eH such that ^ A. 

Proof. We present below 3 concrete examples verifying the above claims. To il- 
lustrate (i) consider the diagonal transformation T4 : M.4[x] R^lx] defined by 
the sequence (Aq, Ai, A2, A3, A4) = 153, 2M' sUs)- determinant cri- 

teria (3) of Theorem 13.51 the operator T4 preserves positivity. However, one can 
check that its inverse sends the non-negative polynomial {x+ 1)^ to the polynomial 
{x + l)(305a;'^ -I- 495x^ -I- 243a; 4- 29) possessing two real and two complex roots. 
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This example shows that in the finite-dimensional case there is a diagonal trans- 
formation which preserves positivity, but whose inverse docs not preserve hyper- 
bolicity. We can extend this example to the infinite-dimensional case as follows. 

By [U Proposition 3.5] that there exists an infinite sequence {Xi}°Za ^ ^ such 

that the sequence of inverses ^ t- f ^ As an explicit example one can take 



P + 5i'^ + 33i + 29' 

An example illustrating (ii) is given in p. 520], see also [3 Example 1.8]. 

Namely, the sequence {1 + i + i^}^o corresponds to a diagonal transformation 

. ( 1 \°° 
preserving hyperbolicity. However, the sequence of inverses < — > leads 



1 + i + i 



=0 



to a diagonal transformation which is not a positivity-preserver. □ 

Definition 3.8. We say that a diagonal transformation T^, a G NU oo, generated 
by the sequence {A^jf^g is a complex zero decreasing sequence (CZDS for short), if 
for any polynomial p{x) G RQ[a;] the polynomial T{p) has no more non-real roots 
(counted with multiplicities) than p. We denote the set of all CZDS by TZ. 

Remark 3.9. Obviously, any CZDS preserves hyperbolicity, i.e., TZ C H. For a 
while it was believed that TZ = Ti. until Craven and Csordas found a counterexample 
[1]. Additionally, one can see directly from the definition that the inverse of any 
positive CZDS is a A-sequence, that is, a diagonal positivity-preserver. 

Finally, Theorem 4.6.13 of [S| states that T e A if and only if T"! € TZ, which 
we disprove below. 

Proposition 3.10. There exist T G A such that T^^ ^ TZ. 

Proof. Use the first two counterexamples from the proof of Proposition 13. 71 □ 

4. Linear ordinary differential operators of finite order 

Our aim in this section is to prove Theorem [XI i.e., to show that there are no 
positivity-, non- negativity-, and ellipticity-preservers which are linear differential 
operators of finite positive order. In fact we are going to show that for any linear 
differential operator U of order A: > 1 there exists an integer n such that U : K„ [x] 
M.n[x] is not a non-negativity preserver. Moreover, we show that one can always 
choose n = 2k. Since any positivity-preserver is automatically a non-negativity- 
preserver and any ellipticity-preserver is a positivity-preserver up to a sign change 
we will get Theorem 1X1 in its complete generality from the above statement. 

Denote by S[si, . . . ,Sk] the ring of symmetric polynomials with real coefficients 
in the variables si, . . . , s^. Let cr; be the Z-th elementary symmetric function, i.e., 

ai = ^ Sj^ ■ ■ ■ Sj^ e S[si, . . . ,Sk], l = l,...,k. 
ji<...<ji 

We will need the following technical fact. 

Proposition 4.1. Let p[x) — {x — xi)^ ■ ■ ■ {x — Xk)^ G R[x, xi, . . . ,Xk]. Consider 
the following two families of rational functions: 

wi = wi{x,xi, . . . ,Xk) = — 1 = 1,. ..,k; 

p{x) 

ui=cri[ — - — — - — y l = l,...,k. (6) 



X X"]^ X k 



Ther 
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1. WieS 



1^1 k. 



2. For any I — 1, . . . , fc one has that 

wi = 2l\ui . . . 

where gi € R[yi, ■ ■ ■ , yi-i], I ^ 1,2, . . . , k, are certain polynomials (that can be found 
explicitly but we will not need their explicit form; in particular, gi = 0). 

Proof. Set pi{x) = {x — xi) ■ ■ ■ {x — Xk)- Then one can immediately check that 

Pi (a;) 

for all i = I, . . . ,k. Using the Leibniz rule we get 

= (p?(x))W = 2p^p{x)p,{x) + ^ c,y;Hx)p['\x), 

id>T-,i+j=l 

where Cij > are certain binomial coefficients. Thus 

i,J>l,i+J=l > I i,j>l,i+j=l 

The result follows. □ 

We are now ready to prove the first main result of this paper. 

Theorem 4.2. Let Uq : R[x] M.[x] be a linear ordinary differential operator of or- 
der k > 1 of the form {IJ) with polynomial coefficients Q ~ {qQ{x), qi{x), . . . , qk{x)), 
Qi G R[x], qk{x) ^ 0. Then for any such coefficient sequence Q the operator Uq 
does not preserve the set of non-negative polynomials of degree 2k. 

Proof. We assume that Uq ^ 0. Since Uq{1) — qo{x) an obvious necessary con- 
dition for the operator Uq to preserve non-negativity is that go (a;) itself is a non- 
negative polynomial. Moreover, qo{x) does not vanish identically by Lemma |2. II 
We will now construct a non-negative polynomial 

p{x) := Pxt....,x^ix) = (a; - xi)'^{x - X2Y ■ ■ ■ {x - Xkf G K[a;,a;i, . . . ,Xk]. 

such that its image under the action of Uq attains negative values. For this we 
define 

R{x,x,,...,Xk) = " = q„{x) + qi{x)i^ + q2{x)i^ + . . . + qk(x)i-^. 

p(x) p{x) p{x) p(x) 

Then in the notation of Proposition 14. II we have that R(x, xi, . . . , Xk) — qa{x) + 
qi{x)wi + . . . + qk{x)wk. Let us fix a;o G K such that xq ^ and for any i = 
0,1,..., /c either qi(x) = or qi{xQ) ^ 0. Set Ui — qi{xo), i — l,...,k, and 
r{xi, . . . , Xk) = R{xq, xi, . . . , Xk). Then r{xi, . . . , Xk) is a linear form in w'l, . . . , w'f., 

k 

where w[ — Wi{xQ). Thus there exist ai, . . . , G M such that qo{xQ) -f X] otiO-i < 0- 

(Notice that by our choice of xq one has go(a^o) > 0.) 
Let now bi, . . . ,bk G M be defined by 

6i=ai, fe; = ^(a— 5i(&i,...,6j_i)), i = 2,...,k, (7) 

where gi are defined in Proposition 14. II Consider the system of equations 

1 1 

XQ~tl 'xQ-tk, 



l,...,k, (8) 
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with unknowns ti,...,tk- It follows from the Vieta theorem that the fc-tuple 
(ti, . . . jtfc) solves system (IHI) if and only if \ , . . . , j_ are the roots of the 
equation 

Note that since &i e M the roots of the latter equation are either real or com- 
plex conjugate. Wlog we can always assume that they are ordered so that zi = 
zi, . . . , 2:21-1 = zi7 G C, Z2i+i, ■ ■■ ,Zk e M. Thus the fc-tuple (ti, . . . , ife), where 
ti = xq — J-, i = 1, . . . ,k, solves system Obviously, ti — t2, ■ ■ ■ , ^21-1 = G C, 
t2i+li ■ ■ ■ ,tk G 

Let us substitute a; = xq, — ti, i — 1, . . . , /c, in the functions wi, I — 1, . . . , k, 
defined in Proposition 14.11 By the definition of the u^'s and using the fact that 
the ti's solve system ^ we get that Ui{xo, ti, . . . ,tk) = h, i = 1, . . . , k. Thus ([7]) 
implies that 

Wi{xo,ti, . . . ,tk) = ai, i = l,...,k. 

Hence R{xo,ti, . . . ,tk) < by the choice of the a^'s. Thus Ua{pti,...,tk){xo) < 
since Pxi,...,Xk i^) ^ for all a; S M and all Xi,i = 1, . . . , k. This contradiction proves 
Theorem 11121 □ 



As we mentioned above Theorem 14.21 together with the results of §2 imply The- 
orem El 



We are now going to strengthen Theorem 14.21 and show that wide subclasses of 
linear ordinary differential operators of finite order k do not preserve non- negativity 
even in degrees much smaller than 2k. In particular, the next statement shows that 
no linear differential operator of odd order k preserves the set of non-negative 
polynomials in 

Proposition 4.3. Let Uq : R[x] — > M.[x] be a linear differential operator of odd or- 
der k > 1 of the form (OP with polynomial coefficients Q = {qo{x), qi{x), . . . , qk{x)), 
qi{x) € ]&[x], i = 1, . . . , A:, qfc(x) ^ 0. Then Uq does not preserve the set of non- 
negative polynomials of degree smaller than or equal to k -\- 1. 

Proof. Consider the polynomial pt{x) — {x — t)'^^^. It is non- negative since k is 
odd. Note that 

rQ{x, t) := ErM^ ^ + (k + l)qi{x)^ + . . . + {k + l)lq,(x)-^-^ 

and set u := ^-7. Then the function 

x — t 

g{x,u) qnix) + {k + l)qi{x)u + . . . + (fc + l)lqk{x)u'' 

is a polynomial in u. Fixing xq such that qk{xo) ^ we obtain that g(xQ,u) is a 
polynomial in u of odd degree. Hence there exists uq such that g(xQ, uq) < 0. Now 
for <o = 2-0 - we get that rQ(a;o,to) < 0. Thus J7Q(pto)(2^o) < since pt{x) > 
for all {x,t). This contradiction finishes the proof. □ 

The next result shows that there is also a large class of linear differential operators 
of even order k which does not preserve non- negativity in M/c[a;]. 

Proposition 4.4. Let Uq : M.[x] M.[x] be a linear differential operator of even 
order k of the form with polynomial coefficients Q = {qo{x), qi{x), . . . , qk{x)), 
qi{x) G W^x], i — 1, . . . ,/c, qk{x) ^ 0. Assume in addition that either there exists 
G R such that qkixo) < or there exists Xq G K such that qk{xa) — and 
qk-i{xo) 7^ 0. Then Uq does not preserve the set of non-negative polynomials of 
degree smaller than or equal to k. 
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Proof. The polynomial pt{x) — {x — t)'' is non-negative since k is even. Similar to 
the above one has 

, , Uniptix)) / \ , / \ 1 , I / \ 1 

As before we set u := ^-r and consider the function 

x — t 

g{x, u) := qoix) + kqi{x)u + . . . + k\qk{x)u^ , 

which is a polynomial in u. If there exists xq such that qk{xo) < then g{xQ, u) is 
a polynomial in u which is negative for sufficiently large values of u. If qk{xo) = 
then g{xo,u) is a polynomial in u of odd degree. In both cases there exists uq 
such that g{xo,uo) < 0. Now for to = xq — we get that rQ{xo,to) < 0. Thus 
UQ{ptg){xo) < since pt{x) > for all {x,t). This contradiction accomplishes the 
proof. □ 

Corollary 4.5. Let Uq : W[x] — > M.[x] be a linear differential operator of order k of 
the form (QJ) with polynomial coefficients Q — (qo^x), qi{x), . . . , qk{x)), qi{x) £ M.[x\, 
i = 1, . . . , k, qk{x) ^ 0. Assume that there exists an even integer i G {2, 4, . . . , fc} 
such that cither there exists G K. such that qi^xo) < or there exists G K 
such that qi{xo) = and qi^i{xQ) ^ 0. Then Uq docs not preserve non- negativity 
in M; [x] for any I > i. 

Proof. If I > i and Uq preserves non-negativity in M/[a;] then Uq preserves non- 
negativity in M.i [x] . The restriction of Uq to the space [x] is given by 

d d^ d^ 

qo{x) + qi{x)— + q2{x)-r^ + ... + qt{x) — 
ax dx^ ax* 

and the result follows from Theorem 14.41 □ 



Remark 4.6. In particular, if i = 2 in CoroUarv 14.51 then there is no degree I, 
I ^ 0, such that Uq : M.i[x] R;[x] preserves non-negativity. Thus for a "generic" 
linear differential operator Uq : K[a;] — > M.[x] with non-constant coefficients there is 
no / such that Uq preserves non-negativity in Mi[a:]. However, the following exam- 
ple shows that there are linear differential operators with non-constant coefficients 
which preserve positivity on Mfc[x] for any even k. 

Example 4.7. Let k be an even positive integer and consider the linear differential 
operator of order k given by 

Ul'^ = 1 + (ax^ + • ^k[x] ^ ^k[x]. 

Then U^''' preserves positivity on IR.fe[a;] if a>0, b >0. 

Indeed, U^'''{p{x)) — p{x) + klak{ax'^ + b) > p{x) > for any positive polynomial 
p{x) = akx'' + ... e Rk[x]. 

5. Linear ordinary differential operators with constant coefficients 

In this section we will prove Theorem [BJ Take a sequence a = (ao, ai, . . . ,ak) 
of real numbers. Denote by Ua the following linear differential operator of order k 

d d^ 
Ua ^ ao + ai— + . . . + Uk (9) 

with constant coefficients. 

By Theorem 14.21 there are no finite order linear differential operators on M.[x\ 
preserving positivity. However, in the case of polynomials of bounded degree, i.e., 
belonging to the finite-dimensional space Mfe[a;], there are such linear differential 
operators, see Example 14.71 
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Theorem [B] follows easily from the next statement of Remak 
[7] which for the sake of completeness we present with its proof. 



and Hurwitz 



Theorem C. For an even integer k — 21 and a sequence of real numbers a — 
{ao, ai, . . . , Q!fe) consider the linear ordinary differential operator (0) with constant 
coefficients. Then the operator Ua preserves positivity (resp., non-negativity) in 
^k[x\ if and only if one of the following two equivalent conditions holds: 

.+aix+ao 



(1) for any positive (resp., non-negative) polynomial p{x) 
one has that 



J7q(p(0)) = floao + fliofi + . . . + klokctk > {resp., > 0); 
(2) the following + 1) x (/ + 1) Hankel matrix 



Wax 
2!a2 



l!ai 
2!a2 
3!a3 



2!a2 
3!a3 
4!a4 



{I 



llai 

- l)!a;+i 
-2y.ai+2 



\llai (l + iy.ai+i {l + 2)\ai+2 ... {2l)la2i J 

is positive definite (resp., positive semi-definite). 

We start with the following observation. 

Lemma 5.1. The operator Ua ■ '^k[x] ^k[x] of the form commutes with shifts 
of the independent variable x. In other words, for any polynomial p{x) S IRfe[a;] set 
q{x) — Ua{p{x)). Then for any a;o G K. we have that q{x — xq) = Ua{p{x — xo)). 

Proof. Take any p{x) = a/a;' + . . . + 
have that (x) = a,fta;'"* + . . . + a^i!. Thus ip^''>){x 
the coefficients of [/„ are constant, the result follows. 



I < k. Then for any positive integer i we 
Xo) — {p{x — xo)y'^K Since 



□ 



Proof of Theorem\^ The equivalence between the conditions (1) and (2) in the 
formulation of Theorem [C] is exactly the same fact as the equivalence between (2) 



and (3) in Theorems 



for Aq — aoj = i^-Oii, i = 1, 



, k. It is proven in the 



required generality in Theorem 11, p. 133 of [17|. What we need is to show that the 
assumption that Ua is a non-negativity-preserver (resp., a positivity- preserver) in 
Mfc[x] is equivalent to condition (1). Indeed, if Ua preserves non-negativity in ]Rfe[a:], 
then ooao -I- aiai -I- ... -I- k\akak — Ua{p(fi)) > for any non-negative polynomial 
p{x) £ Mk[x]. Assume now that for any non- negative polynomial p{x) one has that 
00^0 + aiai + . . . + klakcuk > 0. Set q{x) :— Uq{p{x)). By assumption we have 
that q{0) > and we want to show that q{x) is non-negative. For any xo G M 
consider gxo{x) := q(x -\- xo). By Lemma [01 we have that gxo(x) = Ua{p{x + xq)), 
but /(x) := p(x-|-xo) is a non-negative polynomial. Thus by condition (1) one has 
Uaifm > 0, i.e., 

qixo) = gxM = Ua{fm>0 

for any xq G M. Simple additional consideration shows that the same argument 
with the strict inequality in condition (1) works for positivity-preservers. □ 

Remark 5.2. Theorem [Cl provides the classification of linear differential operators 
with constant coefficients of an even order k which preserve positivity in [x] . On 
the other hand, by Theorems El and 14.21 there are no linear differential operators 
with constant coefficients of even order k that preserve positivity in R2fc [x] . Below 
we bridge this gap between k and 2k for operators with constant coefficients by 
showing that there are no such operators of order k that preserve positivity (or 
non-negativity, or ellipticity) in Mi[x] for any I > k. 
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Proposition 5.3. Let k he a positive integer and let a — (ao, ai, . . . , ctfe) he a 
sequence of real numbers. Consider the operator Ua of the form Then for any 
I > k the operator Ua '■ [x] ~* R/ [x] does not preserve positivity. 

Proof. Wlog we can assume that I is even. We can also assume ao > and at 
least one more entry aj in the sequence (ao, ai, . . . , a^) is non- vanishing. (The 
cases when either ao < or only ao is non-vanishing are trivial.) Take any (not 
necessarily positive!) polynomial p(a:) = akX^' + .-.+aix+ao of degree at most k such 
that ao > and Ua{p{0)) — aoao + aiai + ... + k\akak < 0. Since both ao and aj are 
non-vanishing such a p{x) always exists. Consider now P{x) = Mx^ -\- p{x) where 
M is a large positive constant. By our assumptions one can always choose such a 
large M that P{x) becomes positive. At the same time Ua{P{0)) = Ua{p{Oj) < 0. 
The latter contradicts to the condition (1) of Theorem O implying that Ua does 
not preserve positivity in M; [x] . □ 

Remark 5.4. Finally, notice that in order to get Theorem [Bl from Theorem [Cl it 
suffices to observe that for any positive even integer k the action of the operator 
Ua of infinite order of the form ^ on the space Rfc [x] coincides with the action of 
its truncation 0. 
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